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ABSTRACT 

Multipole expansion of spatial three-point statistics is introduced as a tool for investigating and displaying 
configuration dependence. The novel parametrization renders the relation between bi-spectrum and three-point 
correlation function especially transparent as a set of two-dimensional Hankel transforms. It is expected on the- 
oretical grounds, that three-point statistics can be described accurately with only a few multipoles. In particular, 
\ we show that in the weakly non-linear regime, the multipoles of the reduced bispectrum, Qi, are significant only 

up to quadrupole. Moreover, the non-linear bias in the weakly non-linear regime only affects the monopole order 
of these statistics. As a consequence, a simple, novel set of estimators can be constructed to constrain galaxy 
bias. In addition, the quadrupole to dipole ratio is independent of the bias, thus it becomes a novel diagnostic 
of the underlying theoretical assumptions: weakly non-linear gravity and perturbative local bias. To illustrate 
the use of our approach, we present predictions based on both power law, and CDM models. We show that the 
presently favoured SDSS-WMAP concordance model displays strong "baryon bumps" in the Q;'s. Finally, we 
sketch out three practical techniques estimate these novel quantities: they amount to new, and for the first time 
edge corrected, estimators for the bispectrum. 

' Subject headings: cosmic microwave background — cosmology: theory — methods: statistical 

t> : 

Tj" ■ 1. INTRODUCTION 

Three-point statistics are on track to become main stream tools in astronomy. They have been used by several authors with consid- 
■^-j- ' erable success to constrain the statistical bias between t he distribution of g alaxies and dark m atter ( e.g. lFr\lfl994llJing & Boernerl 
O ■ U998HFrieman & G aztanagal l^alSzapudi et alil2000llScoccimarro et all 1200 111 Verde et al.l. r2002). as well as to constrain primor- 
dial non-Gaussianity of the Cosmic Microwave Background (Komatsu et alll2003l) . 

Nevertheless, the three-point correlation function and its Fourier-transform pair, the bispectrum, are complicated objects. In their 
most general form, they depend on three vectors, i.e. nine variables. Even after translational and rotational symmetries are taken 
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O ' into account, three-point functions depend on the size and shape of a triangle, i.e. still three parameters are needed. Exploration and 
,H \ visualization of the full configuration space becomes a surprizingly daunting task, and most previous studies have been forced to 
■ restrict their investigation to a few hand-picked triangle shapes and sizes, whi ch are not nec essarily representative. 

In the past essentially three distinct parametrizations have been proposed. Peebles ( 1980) used the p arameters r,u,v, where th e 
!> . three sides of the tria ngles are r = ri < T2 < r$, and u ~ T2jr\, and v = (Y3 — T-2)jr\. Alternatively. ISuto & Matsubaral ( fl994l) : 
Szan udi et alJ (2001 ) have used logarithmic bins for the three sides of a triangle. Then shapes can be uniquely described by a pair of 
integers {b\ — 63, 62 — 63) formed from the bin numbers b\ < 62 < 63. Finally, several authors parametrized triangles with two sides 
£^ and the angle between then, r\ , r% , 9 (e.g. lFrvlll994l) . 

Systematic exploration of the full configuration space, and comparison of results are somewhat tedious using any of the above 
parametrizations. Moreover, the connection between real and transform space representations (three-point function and bispectrum) 
is somewhat obscure, whichever description one uses. Hoping to alleviate these shortcomings, we introduce the multipole expansion 
of three-point statistics, which arguably expresses the rotational symmetries in the most natural way. The next section presents the 
definition of the three-point multipoles and the relationship between the multipoles in real and Fourier space. Section 3 applies these 
results in the weakly non-linear regime, bias, and presents predictions for power law and CDM models. The final section contains 
summary and discussion of the results, including measurement techniques for the proposed quantities. 

2. MULTIPOLE EXPANSION OF THREE-POINT STATISTICS 

Since the bispectrum can be parametrized with two lengths, fci, &2, and the angle 9 between them, it is natural to define the 
following multiple expansion 

B{k 1 ,k2,6) = ^2B l (k 1 ,k 2 )Pi{cose)—^, (1) 
1 

where P; (cos 9) is the Z-fh order Legendre-polynomial. Given the two scalar lengths, Bi can be obtained through integration Bi = 
2tt J BPid cos 9. The multipole expansion of the real space three-point correlation function, £f, is defined entirely analogously to 
B l . 

This novel choice of parametrization, while in principle equivalent to previous choices, provides a surprisingly new perspective. 
It facilitates predictions, measurements, and their visualizations since it expresses naturally the rotational symmetry inherent in the 
three-dimensional statistics. To begin with, we show that there is a simple transformation between Bi and £ ; 3 . 
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The relationship of the three-point correlation function £ 3 to the bispectrum B is a Fourier transform: 

£(»•!, ra.rg) = J n? =1 rf 3 fc i S(A; 1 , fc 2 , fc 3 )e 4(fcl ^ +fc2X2+ ' c ^ 3) fe(fc 1 + fc 2 + k 3 ), 



(2) 



(5/3 is the Dirac delta function, and temporarily vector notation is used. This is a six dimensional integral and difficult to treat in 
practice. 

If the general relation is rewritten in spherical variables in terms of the above multipole expansion, the exponential can be expanded 
into spherical harmonics while the addition theorem can be used to expand the Legendre polynomial. The angular dependence then 
can be integrated using the orthogonality of spherical functions yielding the following final result 



(3) 



The formula is analogous to that of the power spectrum, where a formally 3 dimensional integral can be rendered one dimensional 
by = / k 2 dk/(2ir 2 )P(k)jo(kr) in the natural variables. 

3. THEORETICAL PREDICTIONS IN THE WEAKLY NON-LINEAR REGIME 

Perturbation theory predicts ( iFrvl 11984 iGoroff et ail Il986t iBouchet etall Il995t iHivon et all 1 19951) that the bispectrum in the 
weakly non-linear regime is 



1 + ~ M ) P (x) + (Ji + Q Px(x) + ~ (1 - n)P 2 {x) } Pi /., )./>! k 2 ) + pmn. 



(4) 



where x is the cosine of the angle between the two wave vectors, Pi are Legendre polynomials, and P(k) is the linear 
power spectrum. The parameter a depends on the expansion history of the universe; its value is well fitted by 3/7f2 -1 / 140 
( Kamionkowski & Buchalter, 1999). We have written the above formula in terms of Legendre polynomials, to make it explicitly 
clear that the first permutation depends only on terms up to quadrupole. The other two permutations, however, will contribute higher 
I terms. Typically these higher order terms drop quickly, although an infrared divergence of the power spectrum can cause oscillation 
for isoceles (fci = fc 2 ) configurations only. 

In the following we focus on the reduced bispectrum Q = Bj (P\Pi + P2P3 + P3P1), for which the dependence on the amplitude 
of power spectrum is scaled out, thus facilitating interpretation. Analogously to the previous definitions, we introduce the multipole 
expansion of the reduced bispectrum, Qi(ki, /c 2 ). Real space version of this quantity Qi(ri, r 2 ) can be introduced exactly the same 
way, however, in what follows, for Qi will denote transform space quantity. In the literature the notation Q3 is often used for the 
reduced bispectrum: since this paper deals with three-point statistics only, the index on Qi denotes multipoles. 

Figure 1. plots the multipoles of the reduced bispectrum, (21 + l)Qi/4ir for 3 configurations, assuming a power-law power 
spectrum with index n = — 1 for = 1. Multipoles higher then I ~ 2 appear to be negligible with the possible exeption of a 
tiny negative octupole. The right panel displays the configuration dependence; Qi = Q;(fci/fc 2 ). Several values of fc 2 are plotted, 
nevertheless they all coincide: the reduced bispectrum depends only on the ratio of the wave vectors for a power law initial power 
spectrum, as expected. Therefore, in the weakly non-linear regime, most of the information contained in the full configuration space 
can be displayed in three graphs, which show the dependence on the ratio up to quadrupole. 

Most of these results hold for CDM power spectra, with the exeption that Qi will depend slightly on the magnitude of the wave 
vectors as well. Next we present predictions for the monopole, dipole, and quadrupole moment of the reduced bispectrum in CDM 
models. 
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FIG. 1. — Left: The multipole expansion of the reduced bispectrum in a power law n = — 1 model is plotted for three different configurations of fc2/fci: 
isoceles 1/1 (thin dash), 2/1 (thick dash) , 4/1 (solid). Right: the configuration dependence of the reduced bispectrum, monopole (thin dash), dipole (thick dash), 
quadrupole (solid). &2 can be fixed at an arbitrary value due to scale invariance of the initial power spectrum (the actual plot shows three sets of five indistinguishable 
curves calculated numerically for the same values of as for Figure 2.). The x axis shows k\/k2, the shape factor. 
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We define our concordance cosmological model based on WMAP (Spereel et al., l200l and SDSS iTegmark et all 120031 
iPooe & etall 120031) with pai-ameters h = 0.685,0 = 0.309, n b = 0.0228//i 2 ,n = 0.966. Figure 2 displays results for a set of 
fixed &2 = 0.01 — 0.05 as a function of 1 < k±/k2 < 10. As expected, the configuration dependence is not a function of the ratio 
only, but is slowly changing with the magnitude as well. Still, virtually all information can be summarized in one plot. The left 
panel uses a BBKS power spectrum (Barde en et all fl986). with no baryonic oscillations. The configuration dependence is smooth, 
qualitatively similar to th e case of the power law power spectrum. The right panel of uses the same parameters with the more accurate 
EH fit power spectrum (Eisenstein & Hul[l9 98): it features prominent baryonic osciallations. Baryonic oscillations are also evident 
on corresponding plots of the bispectrum Bi as well (not displayed). 
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FIG. 2. — Left: Left: the configuration dependence of the reduced bispectrum, monopole (thin dash), dipole (thick dash), quadrupole (solid) is plotted for a BBKS 
power spectrum (no baryons). The x axis shows k\/k,2, and the set of curves display fixed ki = 0.01, 0.02 . . . , 0.05. Right: the same plot is repeated for EH 
power spectra, featuring baryonic oscillations. 



4. BIAS 

The above theory can be used to construct estimators which constrain the bias in the w eakly non-linear regime. Indeed, in the 
weakly non-linear regime one can expand the biased field as f(S) = bo + b\8 + 62/2<5 2 (e.g. Frv & Gaztanaga, 1993), where 5 is the 
purely gravitational dark matter field. Then the biased reduced bispectrum transforms as q = O lb + bilb 1 (lFrvlll9 94). where the 
lower case denotes the galaxy (measured), and the upper case the dark matter (theory) values. It is clear that bi can only effect the 
monopole term. Thus a simple estimator for the bias can be constructed as 

b _ Qi _ Q2 

11 12 

b 2 = qob 2 - Qob. (5) 

According to the equations, the quadrupole to dipole ratio does not depend on the bias, thus it serves as a novel, useful test of the 
underlying assumptions: a quasi-local perturbative, deterministic bias model and perturbation theory. Figure 4. shows the dipole to 
quadrupole ratio for BBKS, and EH power spectra, respectively. The range of fc's to be used for bias extraction can be determined 
from contrasting the measurements with these predictions. Note that scales where baryon oscillations are prominent are barely 
accessible with present redshift su rveys. On smaller sc ales non-linear evolution is likely to modify these prediction based purely on 
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FIG. 3. — Left: the dipole to quadrupole ratio 3/5Q1/Q2 is plotted for a BBKS theory (thick solid line) for k = 0.01. The right panel shows the same using the 
EH fit for comparison, featuring baryonic oscillations. 
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5. SUMMARY AND DISCUSSIONS 

We have introduce the multipole expansion of three-point statistics, which expresses the inherent rotational symmetry in the most 
natural way. In particular, the relationship of real and transform space become entirely transparent via Hankel transforms. As a 
first application, we explored the meaning of the novel quantities theoretically in the weakly non-linear regime. For the multipole 
expansion of the reduced bispectrum, Qi (k\ , k^j, most of the information can be compressed into three functions, monopole, dipole, 
and quadrupole. From these, we constructed a novel estimator for the weakly non-linear bias. Since the quadrupole to dipole ratio is 
independent of the bias, it serves as a novel diagnostic of the underlying theoretical assumptions. This analogous redshift distortions, 
where the quadrupole to monopole ratio plays a similar role. 

Our repackaging of three-point statistics has yielded novel predictions. In particular, the baryonic oscillations appear to be ex- 
tremely prominent in the presently favoured concordance model, suggesting that our estimators will become useful for main stream 
parameter estimation. 

In previ ous work, a formally simila r parametrization have been used to describe projection of primordial non-Gaussianities of 
the CMB (ISoergel & Goldberg! [l999), and, in an identical result, for projecting the weakly non-linear bispectrum on the sphere 
dVerde et al.ll2002l) . These papers use an intermediate quantity nearly identical to our definition of Bi, thus their results are directly 
appicable to the projection of the multipoles on the sphere. 

The measurement of B{ merits special consideration. The most immediate consequence of our theory is that the 9 dependence of 
£ 3 or Bi should be sampled at the roots of Legendre polynomials of order l max , where l max is the desired ang ular resolut i on. Th en 
Gauss-Legendre integration can yield the £f or B\. Starting out i n real space, one ca n use the algorithms of Moore et alJ ( 1200 ll) to 
measure £ 3 , or, in transform space, one can estimate B directly (Scoccimarro, 2000). Either way, and Equation [5] transforms one 
into the other, which ultimately amounts to a novel edge corrected bispectrum estimator. This is significant improvement over all 
previous techniques, where no edge correction have been used for the bispectrum. 

A theoretically interesting possibility is to expand the estimator of the bispectrum into bipolar spherical harmonics 5k 1 6k 2 dk 1 +k 2 = 
^LMhh ^* 2 / (^ 1 ' ^2X^1 ®Yi 2 )lm fe.g.lVarshalovi ch et all fl988). Since the bipolar spherical harmonics are orthogonal, and {Yi® 
Y/)oo oc Pi, one can show that B\ = l) 1 /\/2l + 1. While the bipolar transformation, especially for low l's, can be calculated 

with fast harmonic transform, the simple minded Gaussian quadrature above is likely to be more practical. Nevertheless, this sheds 
light o n how B\ is related t o the symmetries on direct product of two spheres S 2 ® S 2 . This is the three-point generalization of 
SpICE Szapud iet alJ(l2001l) . eSpICE Szapu dT& et.al.1 ( 1200 3) algorithms for fast, edge corrected estimations of the power spectrum. 

To apply these techniques successfully to data in the framework of high precis ion cosology, it will be necessary to extend our 
calculations with redshift distortions (e.g. iKaiseR Il987t IScoccimarro et ali fl999) taken properl y into account. However, at least 
one (suboptimal) estimate of the real space quantities can be obtained from the transverse power ( Ha milton & Te gmark. 120021) . The 
redshift space bispectrum is parametrized by the five parameters B(k\ t ±, fc2,±, &3,_l, feiji , &2,||)> with _L denoting transverse, and | 
paralell quantities with respect to the line of sight in the distant observer approximation. Then the real space bispectrum can be 
estimated from taking fc lt 11 ~ k 2 ,\\ — 0. 

In the highly non-linear regime it is expected that l max > 2 will be needed f or full description, si nce halo models predict strong 
features in the configuration dependence of the three-point correlation function (Takada lfe Jainll2003l) . Nevertheless, it is likely that 
the multipole expansion will be convenient even in this case. 

Several immediate generalizations to the above ideas are in the works. The ab ove investigations wi ll be applied to the halo model 
to extend the theory for dark matter and galaxies in any regime. According to Takad a"& Jainl (|2003), the halo model integrals are 
7 dimensional, and intractable without approximations. This is no longer the case: with our Equation [3] the integrals reduce to 
three dimensions for each multipole. Generalization of the theory is fairly straightforward for projected quantities, such as angular 
clustering of galaxies, CMB, as well as for vector and tensor correlations, such as CMB polarizations and weak lensing. We will 
be exploring practical implementations of the mentioned measurement methods, with special attention to practical caveats, such as 
binning/regularization and Poisson noise subtraction. Multipole expansion is quite natural for redshift distortions, and we aim to 
generalize o ur calc ulation for the three-point multipoles, preferably without making use of the distant observer approximation as in 
ISzalav et alJ {l998). This and the above generalizations will be presented in subsequent publications. 

It is a pleasure to thank Pablo Fosalba, Jun Pan, and Alex Szalay for stimulating discussions. The author was supported by NASA 
through AISR NAG5-1 1996, and ATP NASA NAG5-12101 as well as by NSF grants AST02-06243 and ITR 1 120201-128440. 
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